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ABSTRACT 

The  subject  of  this  paper  is  the  mechanical  behavior  of  rectangular  plates  subjected  to 
a  combination  of  axial  compression  and  lateral  pressure.  Displacements  and  stresses  are 
obtained  from  a  Fortran  code  based  on  the  von  Karman  plate  equations.  The  effects  of 
various  boundary  conditions,  nonlinearities,  and  imperfections  are  included. 
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INTRODUCTION 


Fatigue  life  predictions  require  knowledge  of  the  stresses  in  a  ship  under  its  operating 
conditions  [Sikora,  Dinsenbacher,  and  Beach  (1983)].  Stiffened  plates  are  a  basic  structural 
component  of  ships  and  submarines.  The  mathematical  equations  governing  the  deformation 
of  thin  elastic  plates  and  methods  for  solution  of  these  equations  are  well-known  [Timoshenko 
and  Woinowsky-Krieger  (1959),  Szilard  (1974),  and  Hughes  (1983)].  The  object  ive  of  our 
work  is  to  use  these  known  mathematical  methods  to  predict  the  stresses  in  typical  ship 
plating. 

Plates  are  subjected  to  axial  tension  and  compression  from  longitudinal  bending  of  a  ship 
due  to  wave  loads.  In  addition,  plates  on  the  bottom  are  subjected  to  lateral  water  pressure.' 
Deck  plating  boundary  conditions  may  be  taken  to  be  simply-supported  ,  whereas  bottom 
plating  boundary  conditions  are  more  closely  approximated  as  clamped.  Deck  plating  may 
be  considered  to  have  an  initial  geometric  imperfection,  whereas  bottom  plating  is  bowed  in 
by  the  water  pressure. 

We  consider  a  rectangular  plate  with  length  a,  width  b,  and  thickness  t.  The  structure 
is  subjected  to  axial  force  F,  and  possibly  uniform  lateral  pressure  p  (force  per  unit  lateral 
area  of  the  plate). 
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PLATE  EQUATIONS 

The  hasir  differential  equations  of  nonlinear  shallow  shell  theory  are 

Dx/"*  W  =  p  +  (h.yy  (wo  +  w},:rx  +^,xr  (^’O  +  w),yy  -2  (^o  +  tv),ry 


(i: 


V  ^  [(^^’O  T  ^)ixy  (^’O  T  ^^')ixx  (^0  T  ^')iyy  ]  ( ^O^tt  ^^’Oiyy  )  (2) 

Here  (x,  y,  ~)  are  Cartesian  coordinates  measured  from  the  center  of  the  plate. 

D  =  Et^/12{\  —  u'^)  is  the  bending  stiffness,  E  is  Young’s  modulus,  and  p  is  Poisson’s  ratio. 

y)  is  the  initial  geometric  imperfection  of  the  plate  midsurface  in  the  z-direction,  while 

v){x,y)  is  the  normal  displacement.  ^'(x,y)  is  the  Airy  stress  function.  Commas  denote 

partial  differentiation  with  respect  to  x  or  y;  e.g,  ; 

4  d‘*w  d^w  d'^w 

Sy  W  -  -\-2W,^xyy  +W,yyyy  -  +  ^Q^2Qy2  ^ 

We  let  u{x,y)  and  v{x,y)  denote  the  displacements  of  the  plate  midsurface  in  the  x  and 

y  directions,  respectively.  These  in-plane  displacements  can  be  related  to  w  and  4>  by  using 

the  strain-displacement  relations 

fx  =  U,:r-\-Wo,xW,:r+lw,l 

Cy  =  U,y+U;o,y«^,y 


c-ary 


=  U,y-i-V,:^-\-Wo,xW,y+Wo,yW,x+W,^W 


>y 


and  the  constitutive  relations 


Et 


r(ex  +  I'^y) 


yv„  = 


Et 


^  1  -t/2 


Et 


{ey  +  I/Cx) 


JV  =  _ f 

2(1 

The  membrane  stress  resultants  (forces  per  unit  length  of  side)  are  related  to  the  stress 
function  by 


IV  =  ^  JV  ~  ^  IV 


xy 


'fxy 
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The  bending  stress  resultants  (moments  per  unit  length  of  side)  are  related  to  the  normal 
displacement  by 

M.X  D(w^xx  ),  1XX  ))  Mxy  -^(1  ^'j^ixy  (3) 

The  bending  stresses  (forces  per  unit  area  of  cross  section)  acting  on  the  top  of  sections 
parallel  to  the  yz  or  xz  planes,  respectively  are 


6Mx 

(^x  —  ^2  ’ 


6My 


We  consider  three  different  possible  sets  of  boundary  conditions; 
Loosely-clamped 


On  X  = 

-b-  • 
X2 

W  = 

W,x  = 

N 

i'ixy 

=  0, 

u  - 

=  constant, 

6  Nx{±%,y)dy  ■■ 
2 

=  F 

On  y  = 

X2  • 

W  = 

W,y  = 

N 

‘■'xy 

=  0, 

V  = 

a 

-  constant,  /Jo 

^y(^)±2^  )d^ 

=  0 

Rigidly-clamped 

On  X  = 

+  - 

:  w  = 

W,x- 

N 

•t  “  xy 

=  0, 

u  = 

h 

-  constant, /^t, 
2 

^x(±j,y)dy  = 

F 

On  y  = 

+  i  . 

;  w  = 

W,y  = 

N 

d xy 

=  V  - 

=  0 

(4) 


(5) 


(6) 


Simply-supported 

b 

On  X  =  +1  :  w  =  w,xx  =  Nxy  =  Q,  u-  constant,  j\  {±%-,y)dy  =  F 

a 

On  y  =  +1  :  w  =  w,yy  =  Nxy  =  0,  u  =  constant,  /J|  Ny{x,+^)dx  =  0 

The  classical  buckling  load  of  a  simply-supported  plate  with  an  integral  aspect  ratio 

y  =  1,2,3, ...  is  given  by  —Fcr,  where 
0 
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SIMPLIFIED  EQUATIONS 


For  normal  displacements  iv  which  are  smaller  than  the  plate  thickness  t,  the  nonlinear 
t('rms  in  (l)-(2)  are  negligible  and  we  can  replace  (l)-(2)  with  the  linear  approximation 


Ny  =  constant  C,  N_ 


0 


D  v'*  =  P  +  Nx{wo  +  +Ny{wo  +  w), 


yy 


(T) 


(8) 


For  the  boundary  conditions  (4)  and  (6)  where  the  lateral  edges  y  =  ±-  are  free  to  expand 

or  contract  Ny  =  0,  whereas  for  the  boundary  conditions  (5)  where  the  lateral  edges  y  =  ±- 

yF  .  .  .  .  .  . 

are  restrained  Ny  =  If  p  0,  equation  (8)  may  be  written  in  the  nondimensionalized' 

form 


V  =  12(1  -  +  — — {wo  +  +  w) 

r  cr  ^ 


_  wEt^  _  X 

where  w  =  — — ,  x  =  —,  etc.  The  nondimensionalized  bending  stresses  are  then 

pb‘*  b 

given  by 


6M2;  1 

pb"^  pb"^  2(1  — 

Note  that  in  the  case  Wq  =  0  the  normal  displacement  w  increases  linearly  with  the  pressure  p, 
whereas  in  the  case  p  =  0  the  normal  displacement  w  increases  linearly  with  the  imperfection 
amplitude. 

We  can  investigate  the  importance  of  nonlinearities  by  substituting  the  linear  solution 
w  to  (8)  into  the  right  side  of  (2).  The  additional  membrane  stresses  SxiSy^Sxy  are  then 
determined  by 

F  F 

Nx  =  —  +  tSx  =  —  +  ^lyy  5  +  c*  =  +  C)  ^xy  —  iSxy  =  — ^,xy 

0  0 


{w,XX+l'W,yy),  etc. 
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■J-  V'*  ^>  =  [(u;o  +  w),ly  -{wo  +  (wo  +  w),yy  ]  -  (wo,ly -Wo, XX  WQ,yy  )  (9) 

If  p  7^  0,  Equation  (9)  may  be  written  in  the  nondimensionalized  form 


where  $  = 


V  ^>  =  [{wo^w),ly-{wQ-Vw),-^^  (u;o  +  rn),37;7]  -  {wo,^y -Wo,^Wo,yy) 


The  nondimensionalized  additional  membrane  stresses  are  then  given  by 


Sxt^  ^ 

pV  ~ 


(10) 


Note  that  the  total  longitudinal  stress  at  the  top  of  a  plate  cross  section  is  —  +  +  Sx-,  etc. 

to  'm 


SOLUTION  METHOD 

These  equations  may  be  solved  numerically  by  superimposing  a  Fourier  series  particular 
solution  plus  Levy-type  homogeneous  solutions  with  coefficients  determined  by  the  prescribed 
boundary  conditions.  This  method  of  solution  has  been  recently  used  by  Kwok,  Kang,  and 
Steele  (1991),  Bird  and  Steele  (1991,  1992),  and  Kang,  Wu,  and  Steele  (1993).  A  thoroughly 
documented  Fortran  code  was  written  to  implement  this  solution. 

Normal  Displacement  and  Bending  Stresses 

The  governing  equations  for  an  isotropic  plate  under  axial  compression  and  lateral  pressure 

load  with  no  imperfection  are  from  (7)  and  (8) 

F  vF 

Nx  =  -r,  Ny  =  0  or  — ,  Nxy  =  0 
b  b 

D  W  =  P  +  Nx  W,XX  +Ny  W,yy 

with  the  clamped  boundary  conditions  from  (4)  or  (5) 

On  X  =  ±1  ;  w  =  w,x  =  0 
On  y  =  ±1  :  w  =  w,y  =  0 

The  solution  to  the  problem  above  is  obtained  in  several  stages: 


7 


1)  l  ii'sl  ;i  solution  i/,'i  for  constant  pressure  load  on  a  clamped  plate  strip  with  Ah-  =  Ny  =  0 
is  (jhlained.  This  particular  solution  that  satisfies 


is  of  th('  following  form 

/  f2y\'^Y 

re,  =  .4  h  -  j  j  (12) 

'l'h('  solution  re,  has  zero  rotation  and  displacement  along  y  =  ±-.  From  (3)  the  bending 
moment  My  is  given  by 


My  Dw,yy 


\6AD 

62 


and  Mx  =  i^My 


p  6^ 

( 'oeflici<'nt  A  in  (12)  is  determined  by  substituting  (12)  in  (11):  ^  ~ 


lie 


Fourier  cosine  series  for  (12)  is  given  by 


re,  =  ^  Wn  cos 


niry 


wlu're 


4  fi  nny  pb^ 

1T„  =  -  /  re,  cos  dy  =  —.  , 

b  Jo  b  SD  \{mTy  (utt)® 


—  1  12  \  .  nn 

—  I  sm 


The  jrarticular  solution  re,  has  zero  rotation  but  nonzero  displacement  along  x  =  +- 

2)  In  order  to  make  the  normal  displacement  zero  along  x  —  need  to  compute  a 

correction  term  that  will  achieve  this  goal.  This  new  term  W2  which  is  the  solution  of  a 
rectangular  plate  with  inplane  resultants  and  Ny  is  computed  in  the  following  way:  A 
solution  of  the  form  W2  =  ^  Xn  (a:)  cos  — ^  is  assumed,  where  Xn  is  a  function  of  x  only 


and  IS  (l( 


'termined  from  the  condition  that  W2  satisfies  the  plate  equation 


D  W  =  Nj:  W,:rx  +Ny 


yy 


(13) 
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with  zero  rotation  and  displacement  equal  to  —W\  along  x  =  ±  — .  From  the  condition  that 
H'l  satisfies  (13),  we  obtain  the  solution 

Xn  (x)  =  An  cosh  V^X  +  Bn  sillll  L^iX  +  COsh  U2X  +  Dn  sinll  J/2X, 


where  u\  and  U'l  are  solutions  of 


D 


2 


-  =  0 


(14) 


The  constants  Bn  and  Dn  can  be  taken  to  be  zero  due  to  the  symmetry  about  the  y-axis,  and 

the  remaining  constants  A„  and  can  be  determined  by  imposing  the  boundary  conditions 

along  X  =  ±7.  Then  Wi  +  1^2  is  a  solution  for  a  plate  under  lateral  pressure  p  and  in-plane 
^  * 

resultants  and  Ny  with  zero  displacement  along  all  the  edges  and  zero  rotation  along 

a  -1  I  ^ 

X  =  ±-,  but  nonzero  rotation  along  y  =  ±-. 


3)  In  order  to  make  the  rotation  zero  along  y  =  ±-  we  first  compute  the  Fourier  coefficients 
of  rotation  on  these  edges  and  then  obtain  a  third  solution  W3  =  W31  +  1032  with  moments 
distributed  along  the  edges  corresponding  to  the  desired  rotations.  A  solution  1031  sinusoidal 
in  the  x’-direction  is; 


_  rmrx  f  cosh  pi  y  cosh  P2y\ 

^31  —  /  .  Cim  cos  I  -  ^  p  p  I 

„r;^d  «  Vcoshyif  cosh  ^2  2/ 

where  the  parameters  pi  and  /i2  are  the  roots  of  the  polynomial 


This  solution  is  obtained  in  a  similar  manner  as  for  solution  W2,  except  that  in  this  case 

TTZTTX 

(Oai  is  assumed  to  be  of  the  form  1031  =  E  u.  (y)  cos  -  and  Ym{y)  is  determined  by 

771  odd  ^ 

imposing  the  zero  displacement  boundary  condition  along  y  =  ±-.  As  defined  above,  u;3i 
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has  zero  displacement  on  all  the  edges  and  zero  moment  My  along  x  =  ib-.  The  rotation 
along  the  edge  y  =  -  is  given  by 


roty  ^  -{w:ix,y)\y^k  =  Cj 

m  odd 

^  ^  C\TTi  COS 


The  rotation  along  the  edge  a:  =  -  is 


m'KX  (  b  b 

cos  -  I  — yi  tanh  yi  -  +  f-i-z  tanh  ^2  r 


a 

m'KX 


D\\r 


(15) 


j  ,  >1  ^  .  mTT  ( cosh  ni?/  cosh  U2v\ 

rotx  =  -(n;3i,.)U=|  =  ^  C^n - sin  —  — — ^ - 

n  odd  «  ^  Vcosh  I  coshy2f/ 

The  expression  in  the  previous  parentheses  is  an  even  function  of  y,  so  it  can  be  expanded 
in  a  Fourier  cosine  series: 


rotx=  5:C..^sin^(E^ 


mK 


12nm  cos 


riKy 


(16) 


\n  odd 


where 


D 


12nm 


A  fi  f  cosh/Ziy  cosh  y2y 
cosh  fj,i  I  cosh  y,2  | 


=  -l 

bJo 


nny 

b  I  cos— dy 


Auk  .  nrr 
sm 


2  2 
ri-M2 


The  moment  at  the  edge  y  =  -  is: 


M, 


f^31iyy  ^^3l-ixx^\y—k  —  ^  ^  C\mD  ^  /ij  -|-  /i2^  COS 


mnx 


Therefore 


Clm  — 


=  X!  ^I'nCOS 


Fu 


mKx 


D(-A  +  /2|) 


(17) 


The  solution  sinusoidal  in  the  y-direction  which  has  zero  displacement  on  the  edges  and  zero 
moment  My  along  y  =  ±-  is  : 
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„  niry  f  cosh  Uix  cosh  i'2X 

=  L  T  , 


where  ui  and  U2  are  roots  of  (14). 

Now  the  rotation  at  the  edge  x  =  —  is 

^  2 


rotx  =  -{W32,x)\x=^  -  X]  ^2n  COS 

n  odd 

/  y  ^2n  COS 


niry  (  ^  u  i. 

-g-  ( -r/i  tanh  —  +  1^2  tanh  —  j 


D 


Tin 


(18) 


The  rotation  at  the  edge  y  =  -  is  given  by: 

nn  .  nn  f  cosh  ui  x  cosh  1/2  x ' 

roty  = -[w32,y)L^k  =  Z^<^2n^sin—  - - - ^ - - - ^ 

'  '  .tf,  h  2  Vcosh^il  cosh  1/25, 


=  EC's, 

Tl  odd 


nrr  .  nrr 


21nm  cos 


rmrx 


a  ) 


(19) 


where 


4m7r  .  mTT 

Dilnm  —  Q  ®  ^  ~o 

z 


/ 

f  m-rr') 
1  2  J 

O' 

[•'i  +  (=f ) 

'0 

The  moment  at  the  edge  x  =  —  is 


M^  =  D  {-W32,xx  -l^W32,yy  )  U=f 


E  C2nD  (-2.^  +  .2I)  COS 
nodd  ^ 


X)  F2n  COS 

n  odd 


mry 

~r 


Therefore 


F2n 

£’(-*'!  +"2) 


(20) 


In  terms  of  the  Fourier  coefficients  of  quantities  rotx,  and  roty,  M^,  and  My,  we  have  the 
following  relationships:  From  equations  (16),  (17),  (18)  and  (20)  we  have: 
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rotx]  =  [D.\J12]  *  [A/y]  +  [DM22]  *  [Mr] 


(21) 


wlici'c 

tn~  .  niTT  1 

L)  i\l  \  —  D\2nm  *  *  Sin  —  *  2  [JT 

{DM22)n,n  =  D22n*j^ - ^  *  (Wentity)^^ 

Siiiiilarly  from  (15),  (17),  (19)  and  (20)  we  have: 

[roty]  =  [DM1 1]  *  [My]  +  [DM21]  *  [Mr]  (22) 

where 

nn  ,  mr  1 

DM21  „,n  =  D2lnm  *  —  *  SIXX  —  *  —-5 - 

b  2  D(^i-/zf) 

(DM11)„„.  =  Dnm  *  ^  *  (Identity )„^ 

U(l.l2  —  Hi  ) 

Ec|uations  (21)-(22)  can  be  written  in  the  matrix  notation: 

DMll  DM12  My  _  roty 
DM21  DM22  Mr  rotx 

or 

[DM]  *  [M]  =  rot 

Here  [DM]  is  the  flexibility  matrix  and  the  inverse  [DM]”'  is  the  stiffness  matrix.  To  take 
advantage  of  the  diagonal  matrices  DMll  and  DM22,  we  write  the  following  relations: 

[Mr]  =  [M]”'  *  ([rotx]  —  [DM21]  *  [DMll]”'  *  [roty]) 

[My]  =  [DMll]”' *([rofy]- [DM12]*  [Mr])  (23) 

[M]  =  [DM22]  -  [DM21]  *  [DMll]”' *  [DM12] 

12 


After  (oinputing  the  vector  solution  to  equations  (23),  we  have  Fourier  coefficients  of  and 
and  from  (17)  and  (20)  we  can  compute  the  Fourier  coefficients  of  the  displacement 
(■orres])onding  to  this  bending  moment.  Then  rui  +  W2  —  103  is  the  solution  which  has  zero 
lotation  on  all  the  edges  as  well  as  zero  displacement. 

Additional  Membrane  Stresses 

The  governing  equation  for  the  additional  membrane  stresses  in  a  plate  with  no  imperfection 
is  from  (9) 


Et 


—  w 


w 


,yy 


with  the  boundary  conditions 


On  J  =  ±f  : 

A  ^,yy 

=  0, 

b 

^,xy 

(±|,  y]  dy  =  0,  u  =  constant 

2 

2 

(25) 

On  y  =  ±1 : 

/Ja  ^,xx 

2 

(x,  ±1)  dx 

=  0, 

f  2  ^ 
^,xy 

(x,±^^dx  =  0,  n  =  constant. 

The  basic  idea  behind  the  solution  method  for  calculating  the  in-plane  Airy  stress  function 
0  is  the  use  of  double  Fourier  expansions  for  both  $  and  the  out  of  plane  displacement  w. 
In  Older  to  solve  the  governing  equation  (24),  we  need  to  calculate  the  Fourier  coefficients 
of  u)  from  the  solution  in  the  previous  subsection.  The  Fourier  expansions  for  w  and  4>  have 


the  following  form: 


•»  =  E  E ( v) (^) 

i even  J  even 

*=EE*.cos(^^)cos(^^) 

I  even  J  even 

With  these  expansions  the  boundary  conditions  (25)  for  $  will  be  satisfied.  The  derivatives 
of  w  are  calculated  by  direct  differentiation  of  (26),  and  the  products  of  derivatives  are 
evaluated  at  each  grid  point.  Then  the  discrete  double  Fourier  cosine  transform  is  used  to 
obtain  the  Fourier  coefficients  of  the  right-hand  side  of  (24),  which  are  stored  in  the  Gaussian 
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rurvature  matrix  A'.  The  Fourier  coefficients  of  <I>  then  follow  from 


0  = 


Etk'ij 


for  i  =  2,4 .  j  =  2, 4 


The  in-plane  displacements  u  and  v  are  obtained  from  the  strain-displacement  relations 
An  appropriate  plane  stress  solution  is  added  to  satisfy  the  boundary  conditions  on  ii  and  v 
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RESULTS 


Normal  Displacement  and  Bending  Stresses 

VVe  first  solved  the  linearized  equations  (7)-(8)  with  no  imperfection  and  the  loosely-clamped 
lu)unda,ry  conditions  (4).  The  shapes  of  the  deflected  midsurfaces  and  values  of  the  nondi- 
mensionalized  longitudinal  and  transverse  bending  stresses  for  plates  of  various  aspect  ratios 
and  axial  load  ratios  are  shown  in  the  figures.  Note  that  for  large  -  the  normal  displacement 
stresses  on  lines  y  =  constant  are  nearly  uniform  at  a  distance  from  the  ends  greater 


a,n  cl 


F 


than  the  width  of  the  plate.  Note  that  as  the  axial  load  ratio  —  is  varied  from  positive  to 

r  cr 

ncygative,  the  waviness  of  the  plate  near  the  ends  increases. 

Maximum  normal  displacements,  and  maximum  and  minimum  values  of  the  longitudinal 
a.nd  ti-ansverse  bending  stresses  are  shown  in  separate  graphs.  The  location  of  the  maxi¬ 
mum  longitudinal  bending  stresses  is  always  on  the  centerline  j/  =  0,  while  the  minimum 

CL 

longitudinal  bending  stresses  are  always  at  the  endpoints  x  =  +-,t/  =  0.  The  maxnnum 

transverse  bending  stresses  are  always  on  the  centerline  y  =  0,  while  the  minimum  transverse 

b  .  .  F 

bending  stresses  are  alwciys  on  the  edgelines  y  =  +-.  Note  that  as  the  axial  load  ratio  — 

is  varied  from  positive  to  negative,  the  maximum  stresses  increase  and  the  minimum  stresses 

usually  decrease.  Also,  note  that  for  most  cases  the  minimum  transverse  bending  stress  is 

the  greatest  stress  in  absolute  value. 

F 

All  values  of  —  were  greater  than  the  buckling  loads,  given  in  the  following  table: 

F  cr 


a/b 

1 

1.2 

1.4 

1.6 

2 

4 

oo 

Buckling 
Loa.d  F/ Fcr 

-2.52 

-2.45 

-2.41 

-2.22 

-1.99 

-1.81 

-1.75 

F 

As  a  check  on  these  results,  the  maximum  normal  displacements  for  the  case  =  0 

cr 


(located  at  the  origin  x  =  y  =  0)  were  checked  against  Roark  and  Young  (1975),  the 

transverse  bending  stress  distribution  on  the  line  *  =  0  was  verified  to  be  nearly  parabolic 

foi-  large  -,  and  the  buckling  loads  were  checked  against  Timoshenko  and  Gere  (1961). 
b 
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Transverse  Bending  Stress  My6/(pb  ),  a/b  =  2,  F/F, 


o 

II 

O 


CO  h'  LO 
LO  T-  T-  OJ 
1-  CD  CO  T- 


(d  P  P  CD 


o  o 
1 


COCMT-;OT-;CVJCO'^_uq 
dcDCD  CDCDCDCDCD 


33 
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Normal  Displacement  wEt^/(pb^),  a/b  =  1,  F/F 
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Normal  Displacement  wEt^/(pb^),  a/b  =  4,  F/F, 


XT  CT) 

CO 

O  K  ■!- 

<y> 

CM  C\J  T-  O 

CO  CO 

O  O  O  ^ 

o  o 

oo  o 

o  o 

(d  CD 

lO 


o 


CO 

in 

CM 

LO 

T~ 

IT) 

O  ID 

o 

CM 

O 

T— 

O 

O 

O 

CD 

o 

CD 

o 

CD 

O 

O 

CD 

CD 

CD 

CD 

42 


Normal  Displacement  wEt  /(pb  ),  a/b  =  16,  F/F, 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  1,  F/F 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  4,  F/F, 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  16,  F/F, 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  1 ,  F/F, 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  2,  F/F 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  4,  F/F, 
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Normal  Displacement  wEt  /(pb  ),  a/b  =  16,  F/F, 
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Maximum  Longitudinal  Bending  Stress  vs  Axial  Load  and  Aspect  Ratio 
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Minimum  Longitudinal  Bending  Stress  vs  Axial  Load  and  Aspect  Ratio 
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Maximum  Transverse  Bending  Stress  vs  Axiai  Load  and  Aspect  Ratio 
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Minimum  Transverse  Bending  Stress  vs  Axiai  Load  and  Aspect  Ratio 
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Maximum  Normal  Displacement  vs  Axial  Load  and  Aspect  Ratio 
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Additional  Membrane  Stresses 

We  next  studied  the  importance  of  nonlinearities  by  solving  equation  (9)  with  no  imperfection 
and  the  clamped  boundary  conditions  (4)  or  (5).  The  graphs  show  the  nondimensionalized 
maximum  and  minimum  additional  membrane  stresses  due  to  bending.  Note  that  the  ad¬ 
ditional  membrane  stresses  are  at  least  an  order  of  magnitude  smaller  than  those  in  the 
preceding  section  from  the  linear  equations.  Note  also  that  the  main  effect  of  pieventing 
lateral  displacement  of  the  edges  i/  =  is  to  increase  the  additional  transverse  membrane 
stress,  whereas  the  additional  longitudinal  membrane  stress  is  unaffected. 

For  these  graphs  the  following  parameter  values  which  are  typical  of  ship  plating  were 
chosen: 

E  =  30,000,000  psi 
u  =  .3 
6  =  36  in 
t  =  .5  in 
p  =  18  psi 

The  additional  membrane  stresses  corresponding  to  other  parameter  values  may  be  obtained 
from  these  same  graphs  with  the  vertical  axis  scaled  in  accordance  with  (10). 
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Maximum  Additional  Transverse  Membrane  Stress  Due  to  Bending  vs 


F/F, 


SR 


Imperfections 


Finally,  we  solved  the  linearized  equations  (7)-(8)  with  no  pressure,  an  imperfection,  and  the 
simply-supported  boundary  conditions  (6).  Let  us  consider  imperfection  shapes  of  the  form 

miTx  nny 

iuq  =  cos  -  cos  — ; —  With  in  =  1,3,5,...  and  n  =  1,3,5,... 

a  0 

or 

rn.TTX  .  mry 

roo  =  sin  -  sm  — - —  with  77i  =  2, 4, 6, ...  and  77  =  2,4,6,... 

a  0 

where  m  a.nd  n  are  the  number  of  half  waves  in  the  longitudinal  and  transverse  directions, 
respectively.  For  these  shapes  we  obtain  the  analytical  solution 


Wo 

W  = - ; - 5 - 

Fcr  ( rnb  n^a  \ 

The  bending  stress  resultants  then  follow  from  (3)  by  differentiation: 


M, 


Fcr  [  mb  n^a  ''  ^ 


cr 

~F 


2a  2mb , 


+  1 


etc. 


F 

Of  course,  to  prevent  buckling  w'e  require  >  —1. 

F zr 


An  alternate  imperfection  shape  proportional  to  the  deflection  of  the  clamped  plate  under 
small  uniform  pressure  without  inplane  load  was  also  allowed  in  the  Fortran  code,  but  for 
brevity  we  do  not  report  on  this  case  here. 
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CONCLUSIONS 


We  have  developed  a  Fortran  code  for  calculating  the  displacements  and  stresses  in  rect¬ 
angular  plates  subjected  to  the  effects  of  combined  loading,  various  boundary  conditions, 
and  imperfections.  The  next  step  is  to  extend  the  code  to  a  three-dimensional  assembly 
of  plates  and  beams  more  representative  of  ship  grillages.  Progress  in  the  analyses  of  such 
structures  has  recently  been  made  by  Danielson  et  al  (1988,  1990,  1993,  1994). 
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